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Abstract 

We consider a model of bilayer lipid membrane with bola-lipids. The bola-lipid is modeled by 
linking tails of the hydrophobic chains in the opposite monolayers within bilayer as a first approx- 
imation. A number of thermodynamical characteristics are calculated analytically and compared 
with the ones of a regular membrane without chains linkage. Pronounced difference between lateral 
pressure profiles at the layers interface for linked and regular bilayer models is found. In the linked 
case, the lateral pressure mid-plane peak disappears, while the free energy per chain increases. We 
have also calculated distribution of the orientaional order parameter of linked chains across the 
bilayer, and found it is in contrast with the usual lipids case. 
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I. INTRODUCTION 



Studying mechanisms of changes in the structure, elastic and thermodynamic properties 
of new artificial membranes is of fundamental interest, as well as is important for construction 
of membrabnes with higher stability under extreme conditions Besides, there are also 
important natural examples of the membranes with such properties, e.g. archaeal lipids 
with distinctive molecular structure that allows the microorganisms to maintain membrane 
integrity in harsh environments 2J. We present below first analytical results describing some 
issential bahaviors of bolalipid membrane that are derived using modification of our previous 
model of lipid bilayer [3]. 

As a first approximation to the bolalipid bilayer membrane, in which lipid tails from the 
opposite monolayers interpenetrate, we consider a model with pairwise linked tails of the 
lipids belonging to the opposite monolayers within a single bilayer, Fig. [TJ Our model being 
rather simplistic nevertheless bears an important property of the bolalipid bilayer in the 
form of constrained meandering freedom of the chains ends in the vicinity of the monolayers 
interface. We found important consequences of this constriction: the entropy of the bilayer 
decreases, the free energy increases, and the lateral pressure Ht{z) and orientational order 
S(z) profiles change drastically. 

The two distinct Ht{z) curves, see Fig. [21 can be understood by comparing orientational 
fluctuations of the hydrocarbon segments of the semi-flexible lipid chains in the linked- and 
nonlinked lipid bilayers. These fluctuations can be characterized by an orientational order 
parameter S(z), see Fig. \Q(z is coordinate measuring depth inside bilayer), calculated using 
our model. The fluctuations reach their maximum at the monolayers interface inside a 
nonlinked bilayer, because the chains ends are free there. Hence, the order parameter drops 
at z — L, see the dashed curve in Fig. HI Simultaneously, a maximum of the entropic lateral 
pressure occurs at z — L, dashed curve in Fig. [2j Distinctly, in the bolalipid "bilayer" 
fluctuations are significantly suppressed at the "monolayers" interface due to restriction of 
orientational freedom of the central segments by their peripheral neighbors. Hence, S(z) 
does not drop at z — L, see solid curve in Fig. |H As a consequence, there is no maximum 
at z = L in the Ht{z) dependence drawn with solid line in Fig. [2] for linked bilayer case 
(quatation marks indicate that in purely bolalipid membrane there is actually no strict 
notion of a bilayer consisting of monolayers sliding upon one another along the interface). 
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We present our analytical results derived in closed form for thermodynamical properties of 
a membrane with linked chains in the weak-disorder limit: i.e. thickness of the hydrophobic 
part of a bilayer is comparable with twice the length of a single monopolar lipid chain. In this 
limit we use more complete version of the energy functional entering the membrane partition 
function than developed earlier jjj: besides the bending energy of a chain conformation, we 
included kinetic energy of the lipid chain. We prove that this makes the path integral 
representation of the free energy of the chains uniquely normalizable. 

Our latteral pressure profile in case of no linkage (see Fig. [2j dashed line) is in good 
agreement with the molecular dynamics simulations [fi[o| see Fig. [3j 

The plan of the article is as follows. In Sec. II we introduce a microscopic model of a 
bolalipid membrane and calculate the membrane free energy using path-integral summation 
over the chains conformations. The inter-chain entropic interactions are treated in the mean- 
field approximation. Several thermodynamic moduli characterizing the bolalipid bilayer are 
derived as well and compared with nonlinked bilayer case. An increment of the free energy 
per chain due to linkage is calculated. In Sees. Ill and IV we calculate analytically the 
lateral pressure distribution (profile) across the hydrophobic core of the lipid bilayer and 
make comparison between the cases with- and without linkage. Also calculated is chain order 
parameter that characterizes correlations between the orientations of the chain segments and 
clearly demonstrates an increase of the orientation order in the bolalipid case as compared 
with nonlinked bilayer. 

II. MICROSCOPIC MODEL OF BOLALIPID BILAYER 

The bolalipid bilayer membrane is modeled by linking pairwise the tails of the chains 
belonging to the opposite monolayers. Hence, a couple of linked chains is substituted by a 
single semi-flexible string of length 2L, where L is the monolayer thickness, see Fig. [51 
Correspondingly, conformations of the string as a transmembrane object obey combined 
boundary conditions at the opposite head group regions of bilayer with coordinates z = 
and z = 2L respectively, as is described below in detail. 

By doing so, we take into account the decrease of chain's freedom in the midplane in case 
of linkage in comparison with regular membrane. With bending (flexural) rigidity Kf, and 
with the mean-field approximation accounting for entropic repulsion between neighboring 
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couples of pairwise linked chains (see Fig. EJ), the energy functional of a single string, E t , 
has the form: 



2L 



P K (z) 



K f f d'K(z) 
2 1 dz 2 



B _ 
— R z z 




(1) 



Here harmonic potential U e ff = BR 2 /2, with self-consistently defined rigidity B, describes 
entropic repulsion between the strings, z is coordinate along the string axis, and R(z) 
is vector in the {x, y} plane characterizing deviation of the string from the straight line, 
R 2 = R 2 X + R 2 (see E]). The choice of harmonic potential is justified since we assume 
finite softness of the effective "cage" created by the neighboring lipid chains in the limit 



of small chain deviations. A harmonic potential was considered in earlier work 



111 ] for a 



semi-flexible polymer confined along its axis. The elastic energy treatment of the molecular 
chains, especially well known in the theory of polymers, see e.g. jl2,[l^,Q]- The first term in 
Eq. ([T]) represents kinetic energy of the string, p is linear density of mass: p = m(CH2)N/L, 
where m(Ci?2) is a hydrocarbon group mass, N is the number of hydrocarbon groups per 
chain (for numerical estimates we took iV = 18, see j^]). 

The bending energy term in Eq. ([1]) represents the energy of the chain trans or gauche 
conformations. It contains the second derivative over the z coordinate rather than over the 
contour length of the chain. This approximation is valid provided that deviations from the 
z axis are small with respect to the chain length L: 



2L ~ \L 2 P eff ) 



: 1. (2) 

This limit is opposite to the one considered in the long polymer theory where the second 
derivative in Eq. (CQ) is substituted by the first derivative in the flexible chain approximation. 

Using the functional Eq. (CQ) the chain partition function is found as a path integral over 
all string conformations: 



Z 



exp 



exp 



E((R{z),IL(z)) 
e({R x {z),R x {z)) 



D R X D R X D RyD Ry 



DR r DR 7 



(3) 
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The second equality in Eq. ([3]) holds when the membrane is laterally isotropic and x and y 
deviations can be considered independently. 

To calculate the path integral Eq. ([3]) we rewrite the energy functional Eq. (jTJ using the 
self-adjoint operator H: 

Et = 2 / ^ + Ri(z)HRi(z)dz) , (4) 

i=x,y 

The operator H is obtained after integrating by parts the expression Eq. (OQ) under the 
following boundary conditions for the string that models two linked chains belonging to the 
opposite monolayers (the z-coordinate spans from one head group at z = to another at 
z = 2L). The chain angle is fixed in the head group region: 

R'(0) = 0; R'(2L) = (6) 

No total force is applied upon chain at the head group: 

R'"(0) = 0; R"'{2L) = (7) 

These boundary conditions, as well as the energy functional in Eq. ([!]) differ from the ones 



used to describe a single monolayer of a nonlinked lipid bilayer (compare 



pR~(z) K f fd' 2 R(z)\ 2 B 



where E™ is the energy functional of a single monolayer, and the motions of the chains in the 
opposite monolayers forming a noniterdigitated bilayer are independent. The total energy 
of a bilayer in this approximation is then twice the energy of a single monolayer: 2 x E™. 
We impose the following boundary conditions for a monolayer: the chain angle is fixed, and 
no total force is applied to the chain at the head group: 

R'(0) = 0; R'"(0) = (9) 

No total force and no torque is applied at the free chain end (i.e. at the monolayers interface 
inside the bilayer): 

R"'(L) = 0; R"(L) = (10) 



Finally, in both cases, the free energy of a bilayer equals F = — ksTlii^Z), where Z is 
partition function of a bilayer. Using expressions Eq. (JT]) or Eq. ([8]) for linked or nonlinked 
bilayer respectively, we differentiate the free energy and obtain the self-consistency equation 
in the form (expressed below for linkage case): 

§ = 2L(fl 2 > (11) 

As in (j^, we take into account that hydrocarbon chains of lipid molecules are bulky objects 
that possess finite thickness and introduce an incompressible area of the chain cross section 
A (see Fig. E]). The area occupied by a lipid chain in the bilayer is related to the string 
mean square deviation (R 2 ) by the following formula 
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cL4 = 7r(R 2 )= (VA-v^To , (12) 



where 5 A is the area swept by the string formed with the centers of the chain cross sections. 
In the text below we imply by chain deviations those of a string described by the R vector. 
The self-consistency equation Eq. fill I) combined with formuli Eq. (1121) permits us to find 
the A dependence of the coefficient of entropic repulsion B and finally derive the membrane 
equation of state in a form of pressure-area isotherm (see Appendix A for details). 

To make numerical estimates based on our model of a lipid bilayer we use the following 
parameters values: monolayer thickness L = 15A, chain incompressible area A = 20A 2 , 



To = 300-fT as reference temperature. The chain flexural rigidity is defined as 15] Kj = EI, 
where E m O.QGPa is the chain Youngs modulus [l?J and I = A^/Air is the (geometric) 
moment of inertia. The flexural rigidity can also be evaluated from polymer theory Q 
Kj = ksTlp, where l p « L/3 is the chain persistence length [17[ and ks is the Boltzmann 
constant. Both estimates give approximately Kj ps k B TL/3 at chosen L and at T = T . 

III. LINKED-CHAINS BILAYER: THE FREE ENERGY INCREMENT 

The eigenvalues and eigenfunctions of the operator H defined in Eq. (J5j) obey the following 
equation: 

HR n = Kf + BR n = E n R n (13) 
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Solving this equation with the boundary conditions Eq. ©-([Zj) one obtains: 



E n = B + -^,k n = irn/2,n>l; E = B, (14) 



Rn(z) = c n cos(k n z/L),n > 1; R (z) = J —, (15) 



where c n = \/l/L and A n = 2-nL/k n is the wavelength. Several eigenf unctions are shown in 
Fig. El 

Then an arbitrary string conformation, described with the deviation from the z-axis, 
R x (z, t), as well as its energy are expanded over eigenfunctions R n and eigenvalues E n found 
from Eq. 

R x (z,t) = J2Cn(t)R n\ z )'i 

1 n= ° (16) 

n=0 n=0 

The bilayer partition function is then found as the integral over the coefficients of expan- 
sion C n and conjugated momenta p n = pC n in Eq. ( |T6l) 



7 f°TT ( Pn C 2 n E n \d Pn dC n 
Z x = | exp 



n=0 



2pk B T 2k B T J 2nh 

7 rji I 7 rp (17) 

n k B T p _ -p-r k B T 
h V E„~ I*- Hw* 



where u n = \J E r J p. It is important that the latter expression for uj n in the limit of a 
free string, B = 0, leads to the well known bending waves spectrum of Euler beam 15| : 



ui n = yEIk^/p (with k n = k n /L), as it follows from Eq. (fl4j) and expression for the 
bending rigidity Kf = EI mentioned above. Hence, by including kinetic energy of the chain 
into the energy functional E t we obtain correct dimensionless expression for partition sum 
in the Eq. (117)1 . 

Using Eq. (1171) . Eq. ([3D and F = — k B T In Z we find the following expression for the free 
energy of a bilayer with bolalipids in our model: 

n=0 " 



Using then relation Si nt = 



—dFint/dT we find the following expression for the entropy Si nt : 



/ f) Z?\ n max 



Huj n J \ u n \ dT J v T 



(19) 



Both expressions are valid provided the motion of the lipid chains at room temperature 
T is classical (not quantum), i.e. : ^Vw^Ab^ « 1. The upper cutoff n max in the 
sums is defined by condition that the shortest half- wavelength 0.5\ nmax = irL/k nmaa , of the 
eigenfunction R nmax is not shorter than the Cif 2- m onomer length of the "chain segment". 
Hence in case with linkage n max = 11 and hjJn max /kBT = 0.29. 

The free energy and entropy of the nonlinked bilayer, F non and S non respectively, are 
obtained using the same relations as in Eqs. f|T8|) and (|19p . but with the corresponding change 
of the frequencies spectrum u n that results from the noniterdigitated bilayer conditions 
expressed in Eqs. (I9l-( fl0j) . Using the above relations we calculated linkage-related free 
energy and entropy "cost" as the differences of the respective bilayer free energies and 
entropies in the linked and nonlinked cases. Our results are represented in Fig. [8] and Fig. [91 
In Fig. [8] the free energy increment (per chain) of the order of 5ksT in the bolalipid bilayer 
with respect to the nonlinked one is caused by the corresponding decrease of the entropy 
~ 5k b (per chain), see Fig. Location of the entropy decrease in the linked bilayer can be 
found by exploring the chain's orientational order parameter S(z) defined as: 

S(z) = \ {3(cos*6(z)} - 1) , (20) 

where 6(z) gives distribution of the tangent angle of the chain across the bilayer. Straight 
(ordered) chain possesses 8 = and S(z) = 1. In the limit of small deviations from the 
straight line 9 < 1 considered in our model the order parameter can be expressed using the 
following relations: 

(cos'e(z)} « 1 - (tg*9(z)) = ((R'(z)f) = M £ , (21) 

n=0 n 

so that finally we obtain: 

ffW Ml _ ?^ E (22) 



n=0 n 

Calculated order parameter distributions across the bilayer, S(z), in our model with and 
without linkage are represented in Fig. HI The solid line corresponds to linked chains (mod- 
eling bolalipids), and dashed line is calculated for nonlinked case. It is obvious from the 
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Fig. H] that main difference occurs at the monolayers interface (z = L) inside the bilayer. 
Free chain ends acquire maximal disorder in this region, while linked tails remain quite or- 
dered. Another manifestation of this mid-bilayer ordering phenomenon will be seen in the 
next section in the calculated behavior of the lateral pressure profile inside bilayer. 



IV. LATERAL PRESSURE PROFILE AND PRESSURE- AREA ISOTHERMS 
FOR BOLALIPID BILAYER 

The equation of state of the lipid chains in the bilayer can be derived as follows: 

'BF, 



8A) T < < 23) 
where P t is the total lateral pressure (here has dimensionality of a tension), produced by 
linked hydrocarbon chains. Substituting expression for the free energy from Eq. (fl~8]) into 
Eq. ([23]) one finds: 

n=0 v / T n 

We may consider P t as an integral of the lateral pressure distribution (profile) function, 
lit (2), over the hydrophobic thickness of the bilayer: 



P t = J U t (z)dz. (25) 

In order to find out lit (2) defined this way, it is possible to use the following formal trick. 
Namely, the dependence on aria A of E n arises via dependence of the "potential" B(A), 
that enters operator H in Eq. ( 1T31) . One may in addition formally consider B(A) as being 
^-dependent function. Then, a well known relation from the perturbation theory [3] leads 
to the following equation: 

dE n \ f {8E n \ fdB\ dz _ [ ri2 , ,fdB\ 



saj t j y5B) T yM) T T = J R ^\M) T iz < (26) 

where 1 means unit length. Now, substituting Eq. (j26|) into Eq. (1241) we find analytical 
expression for the lateral pressure profile from the relation: 

fl= -/ w Sf © f is /«'> & (27) 



9 



Hence, finally: 



Calculated in our model lateral pressure profiles for the bilayer with and without linkage are 
presented in Fig. §2§. It is remarkable, that lateral pressure peak at the nonlinked monolayers 
interface, as seen in the dashed curve, disappears in the linked chains case. Hence, entropic 
repulsion between the lipid chains is indeed weaker in the region where the entropy related 
with the chain orientation order is smaller (compare with Fig. H]). 

Next, it is straightforward to check that due to ortho normality of the eigenf unctions 
Rn{z) the integral of Ht{z) over dz across the bilayer thickness leads again to the expression 
in Eq. (124j) for the total lateral tension P t : 

where we used relation that follows from Eq. ( fT4l) : 

In Figure (TTUT) the calculated pressure-area isotherms for bolalipid (linked chains, solid line) 
and nonlinked (dashed line) bilayer are presented. It is obvious from the figure that lateral 
entropic repulsion responsible for the lateral pressure in the hydrophobic part of the bilayer 
is weaker in the bolalipid bilayer comparatively with nonlinked bilayer at the one and the 
same area per lipid chain and other parameters fixed. 

Differentiation of Pt{A) gives the area compressibility modulus 

(31) 

as a function of the area per chain and temperature. The equilibrium condition is found 
by equating the pressure produced by linked chains to the effective lateral pressure in the 
bilayer: 

P t {A{T)) = P eff = j + P HG + P vdW , (32) 
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where 7 is the surface tension at the hydrophobic-hydrophilic interface; Pug is the head 
group repulsion of electrostatic origin; P v dw is the pressure arising from the vander Waals 
interactions between chains, etc. We choose P e ff > 7 ~ 70 dyn/cm because attractive 
dispersion interactions between hydrocarbon chains are included in the effective surface 



tension 



6|. At room temperature 
one has: 50 < P e ff < 150 dyn/cm 
of linked chains (bolalipids) : 



or a typical lipid bilayer with effective surface tension 
Analytical solution for the total presure in case 



P > M = 3^W- 1)5 /3 • - ^ + >) (33) 

Analytical solution for the total pressure in case of nonlinked chains : 

Oh T 

It follows from the analysis of these expressions that the linkage effect on the total lateral 
pressure at a given area is more pronounced at larger areas per lipid (lower pressures) region, 
that corresponds to hihger orientational disorder of the chains. 

Coeffecient of area expansion Eq. ( 13TT) calculated using Eq. ( 133]) for bolalipid membrane 
is smaller than obtained from Eq. (134"1) for membrane without tails linkage, the latter is in 

n 

close agreement with the measured ones [201 ] . 

In Fig. [11] the temperature dependence of the area per chain in the bilayer is shown. This 
curve is increasing with temperature due to a more frequent collisions of chains. Fig. [12] 
displays Pt{A) dependence. This curve is decreasing with the area per chain due to the 
following reason: when the chains occupy more space they collide less frequently and produce 
less entropic pressure. 

Now we can verify the exploited approximation of the small chain deviations in the 
bilayer Eq. ([2]) . We calculate the thermodynamic average of the chain fluctuation amplitude 
(R 2 (z)) using the relation (i?^ (z)) = J2 n ^Pn) ^n( z ) an d averaging over C n : 

(R 2 ( Z ))=k B Tj2^ 1 - (35) 

n n 

It is worth mentioning that integration of both sides of Eq. ( 135]) over z from to 2L 
provides the self-consistency equation Eq. ( ITT]) . Since E n oc n 4 , the sum in Eq. ( 135]) converges 
fast and allowing for the relation P^{z) ~ 1/L, we can estimate it as J2 n ^/E n oc 1/B. 



11 



According to Eq. (JTj), the increase of potential energy associated with the increase of area 
swept by the string from to 5 A is of order B2L5A (the string is formed by the centers of 
the chain cross-sections). On the other hand, it is equal to the work against the pressure 
P e ff needed to increase the area per couple of linked chains in the bilayer from A to 
A: B2L5A pa P e ff(A — A ). From the last equality and relation Eq. (j!2p it follows that 
B > P e ff/L. Then we evaluate: 

Ai "e/f 
n J J 

and find a rough estimate for the upper limit of the small parameter: 



VWR>/2£ < {k B T/L 2 P ef} )^ = 0.16. (37) 

Finally, we compare the amplitudes of linked chains fluctuations in the bilayer Eq. (133]) 
and in empty space. For a free couple of chains with flexural rigidity Kf the characteristic 
deviation R° can be evaluated by equating the chain bending energy to k B T . This yields 

R0 K kBT L 3 „ L 2 (3g) 

K f 

Allowing for Eq. §E§ and Eq. ([38]), we find ^/(R 2 ) /R -0.1. 
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APPENDIX A: SOLUTION OF SELF-CONSISTENCY EQUATION 



Here we present the solution of the self-consistency Eq. ( TTTj) and find the analytical 
temperature and area per couple of chains dependence of the lateral pressure produced by 
the linked chains using the equation of state Eq. (I23p . It is convenient to perform the 
derivations in dimensionless parameters, 

a = A/A ,b = ^B, (Al) 
K f 

and to introduce the auxiliary parameters 

k n = {Ttn/2)\n>l-v = -^-, (A2) 

where L ~ lb A is the chain length, A ~ A 2 is the "incompressible area" of the chain 
cross section, and the chain flexural rigidity Kf = IcbTL/3 at T « T = 300 K. Using these 
estimates we obtain v = 0.009. 

In the introduced notations Eq. (1A1|) and Eq. flA2[) with E n defined in Eq. (1140 the 
self-consistency equation Eq. (fTTT) acquires the form 

2v(Va-l) 2 . (A3) 



b ^ ^ b + hi 

n=l 11 

The terms in the sum on the left hand side of Eq. (1A3I) decrease fast with growing n and 
we can use integration instead of summation over n. For example, for the effective tension 
P e ff = 70 dyn/cm, we have b ~ 10 3 , while kf ~ 6. In this regime we can solve Eq. (1A3I) 
analytically by substituting summation over n with integration, which yields 

1 f~ in 1 



where we took integral using complex functions theory, and k n is defined in Eq. flA2j) . 
In case of membrane with no linkage (see Fig. \5§ Eq. (1A3|) takes the form: 



n=l 

'e// 

membrane with no linkage Eq. (1A4|) takes the form 



Since b ~ 10 3 {P e ff = 70 dyn/cm), we (approximately) integrate over n, so that in case of 
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V 1 ^ 1 r dn _ 1 



where L, is defined as in 



3j: k n = Tin — 7r/4. 

In both, Eq. (1A4|) and Eq. (IA6I) . we omit 1/6 as 6 ~ 10 3 , which leads to the same 6(a) 
dependences in both cases: 



(A7) 



4w 4 / 3 (v^- 1)8/3 ' 

which is then used in the equation of state Eq. fl23l) . As a result we find the expression 
for the lateral pressure produced by the linked hydrocarbon chains Eq. (|23|) 
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FIG. 1: a - Membrane with bolalipids. b - Our model of membrane with bolalipids. 




FIG. 2: Lateral pressure distribution in the hydrophobic core of the bilayer with (solid line) and 
without linkage, z is coordinate along the chain axis normalized by the monolayer thickness L and 
spanning from one head group [z = 0) to another [z = 2L). The parameters for the lipid bilayer 
are as follows: monolayer thickness L = 15A, area per chain Aq = 20A, chain flexural rigigity 
K = IvbTL/3, temperature T = 300K, pressure normalized by total bilayer pressure P t ° = 140 
dyn/cm. 



1(3 




-20 20 

distance from bilayer center (A) 



4CJ 



7[ correspondingly. 



FIG. 3: Figures a) and b) demonstrate molecular dynamics results from f| and 
In both figures the pressure profile related to the hydrophobic region of membrane is bounded with 
dashed rectangular. It is in a good qualitative correspondence with our results (see Fig. [21 dashed 
line) . 
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0,0 0,5 1,0 1,5 2,0 



z/L 

FIG. 4: Order parameter in cases with bolalipids (see Fig. [1]) and no linkage (see Fig. [5]). 




bi layer 



FIG. 5: Model of lipid membrane in the mean-field approximation: we substitude interaction 
between neightboring chains by an effective quadratic potential. 
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FIG. 6: Hydrocarbon chain as a flexible string of finite thickness. R(.z) is the vector characterizing 
the deviation of the center of the chain cross section from the z axis, |R(z)| = ^ R^(z) + Ry(z); Aq 
is the incompressible area of the chain cross section; A = ir (R- 2 ) is the area swept by the centers 
of chain cross sections; A is the average area per lipid chain in the bilayer. 
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FIG. 7: The eigenfunctions R n {z) of the self-adjoint operator H for the boundary conditions Eq. (J6j) 
and Eq. (|7|). Other parametrs are as in Fig. [2j 
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FIG. 8: Not normalized free energy (per chain) difference of membranes with linkage and with 
no linkage between the opposite chains ends. Membrane with linkage "costs" more free energy at 
fixed area per chain. 
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FIG. 9: Entropy (per chain) difference between membranes with linkage and without linkage. 
Entropy of membrane with linkage is lower than that of a membrane without linkage. 
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FIG. 10: Total lateral pressure comparison. With the same area per chain, membrane with linked 
chains produce less pressure than nonlinked membrane. Area and pressure normalized by constants 
described in Fig. [2] caption. 




FIG. 11: Tempereture dependence of equilibrium area per bolalipid chain, A. Temperature is 
normalized by Tq = 300-fT, area is normalized by Aq = 20 A 2 . 
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FIG. 12: Calculated total lateral pressure Pt produced by linked hydrocarbon chains as a function of 
area per chain at two temperatures T\ (solid line)< T<i (dashed line). Area and pressure normalized 
by constants described in Fig. [2] caption. 
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